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We point out that the free single-fermion propagator which is used in the QFT equations for 
two-fermion systems, has a bosonic structure, transforms to the single-boson propagator for the 
Klein-Gordon equation in the nonrelativistic limit, and, therefore, does not satisfy the Lehmann 
spectral representation for fermions. Proceeding from the Lippmann-Schwinger integral equation, we 
obtain a two-fermion relativistic wave equation in which the free two-fermion equal-time propagator 
satisfying the Lehmann representation, is used. This wave equation has been applied to the electron- 
positron system. In addition to the eigenstates for positronium, we find a solution which describes a 
deep bound state with the binding energy ~ 2m and anomalously small mass of the bound system. 



I. INTRODUCTION 



Calculations of bound states of two-body systems have been approached by several methods. In the Hamiltonian 
approaches the fermion-antifermion or two-fermion wave function satisfies two wave equations or 4 x 4 matrix Dirac 
^ ■ equations which are generalizations of the Dirac equation relative to each particle of the interacting system (see [1-3] 
^sO \ and references therein). 

Systematic descriptions of the two-body system are based on quantum field theory (QFT). Here the starting point 
is the Bethe-Salpeter or Dyson-Schwinger equation (see [4-14] and references therein). To solve the equation analytic 
methods have been used by reducing it, with appropriate hypothesis, to a three-dimensional wave equation. In this 
respect the quasipotential approaches involve the use of delta- function constraints on the relative four-momentum [7] . 
^^ , It is well known that the obtained results for bound states of the electron-proton and electron-positron systems 

""^ ■ are close to the spinless case [6,15] and reproduce the eigenvalues of the Schrodinger equation with the Coulomb 
interaction potential which are corrected by certain relativistic and spin contributions ~ a (a = f- = jkj is the fine 
structure constant) [16,17]. 

Although the eigenvalues obtained for the electron-positron system, are close to the spinless case [6,15], the ap- 



& ■ proaches reproduce two important effects, namely, that the singlet state is more bound that the triplet, and the 
fine-structure splitting has the correct order of magnitude a 4 (a — j- = -^ is the fine structure constant) [16,17]. 
Here we note the following circumstance. The Bethe-Salpeter equation contains the free two-fermion propagator 
. ,-h \ which is reduced to the product of the two free single-fermion propagators, as shown, for example, in [7]). In QFT 
S^ ■ the free single-fermion propagator is defined as: 

H ' 
. C3 ; G°(x - x') =-%< 0|TV>(2#V)|0 >, 

where |0 > is the vacuum ground state; tp(x) is the fermionic field operator and ip c (x) is the conjugate operator. 

At present, the Dirac conjugation is commonly used in QFT: ip c (x) — i^*(x)^°. Then this free single-fermion 
propagator coincides with the Green function for the Dirac differential equation [18]. The latter is defined as: 

( 1 p-m)G°{x-x') = 8 {A) {x-x'), 

where 7 denotes the Dirac matrixes in the standard representation [16] and m is the mass of the particle. 

Thus, the free single-fermion propagator in QFT is adopted to be identical with the Green function for the Dirac 
equation. This propagator is given by: 

G°( P ;E) = ^±^G° KG ( P ;E), (1) 

2m 

where G° KG (p; E) is the Green function for the Klein-Gordon differential equation which describes bosons: 

G° KG (p;E) = -=L (^^y - ET^pj) = Ei-lipf (2) 



Here e(p) = (m 2 + p 2 ) 1 / 2 . 

Note the minus between the two pole functions in Eq. (2). It is due to the commutative relations for bosons. The 
same structure of the Green functions exists for all Bose particles, for example, for phonons in solids, of course with 
account of differences in pole amplitudes and, respectively, in spectral densities. 

In principle, the free fermionic propagator (1) has also the boson-type structure in the sense of the minus between 
the two pole functions. In the nonrclativistic limit when p — > and e(p) — ► to, we have: 

7P + m = 2m f Q Q 

and the free singlc-fermion propagator (1) transforms to the free single-boson propagator (2). It is not unexpected 
because any solution for the Dirac equation is also the solution for the Klein-Gordon equation. 

Thus, there is a very strange result in QFT that the free particle propagator does not depend on particle's statistics 
in the nonrelativistic limit. 

However, according to the Lchmann spectral representation [19], the pole amplitudes of the fcrmion propagators 
must be positive in all regions of the pole spectrum. Unlike boson propagators, the exact singlc-fermion propagator 



can be presented in the form ([19] and see also [20]) 

f°° 
G(p; E)= duj 

Jo 



A(p;u) B(p;w) 



E - u) + i0+ E + w-iO^ 



where A and B are real and positive. 

Thus, we point out that the bare single- fermion propagator used in the QFT approach to the two- fcrmion problem, 1) 
has the bosonic structure; 2) transforms to the boson propagator for the Klein-Gordon equation in the nonrelativistic 
limit; 3) does not satisfy Lehmann's representation for fermions. 

In the present paper we attempt to study the bound states of the two-fermion system. Proceeding from the 
Lippmann-Schwinger integral equation, we obtain a two-fermion rclativistic wave equation in which the free two- 
fermion equal-time propagator satisfying the Lchmann representation, is used. We apply this theoretical formalism to 
the electron-positron system. We omit spin interactions in the system which give small rclativistic and annihilation 
corrections ~ a 4 to the energy levels of positronium. The main result, we would like to point out here, is the possibility 
of unusual bound eigenstates for the electron-positron system. We show that in addition to the positronium states, the 
wave equation describes the deep bound state with the binding energy ~ 2.00794 x m and the exponential damping 
radius of the wave function 1.37812 x A e (A e is the Compton length of the electron). It is essential that similar deep 
bound states cannot exist for two interacting spinlcss particles. Units c = H = 1 are used. 

II. TWO-FERMION EQUATION FOR BOUND EIGENSTATES 

We consider a system of interacting two fermions with the masses m,i and charges g,, i = 1,2. Let (f)(r\,r 2 ;E) be 
the wave function of a bound state of the two-fermion system. Then, according to quantum theory of scattering, this 
wave function should satisfy the Lippmann-Schwinger integral equation [21]: 

0(r 1 ,r 2 ;ii;)-ydr' 1 ^G o (r 1 -r' 1 ;r 2 -r' 2 ; J B)U(r' 1 -r^^(r' 1 ,r^;ii;) (3) 

where V is the interaction between the particles; E is the total energy of the system; G° is the Fourier component of 
the free two-particle equal-time propagator: 

G°(n - ri; r 2 - r' 2 ; E) = J d(t - t')G (x u x 2 ;x' 2 ,x' 1 )\ ti=t , tl - t , ti=h2 , (4) 

where 

G (xi,x 2 ;x' 2 ,x[) = -i< 0|TV'(i)(xi)V'(2)(a;2)V'p)(a;3)V'( l i)(a;4)|0 > . (5) 

Here ip^(x) is the fermionic field operator for the i— particle which can be written as [16]: 

^wW = ^2[(H,p<t>i,p( x ) + b t P ^p( x n- ( 6 ) 



For free wave functions we use: 



<t>i, v {x) = \/— 7 2 Texp(*pr-«e»(p)t)u 1 Tj 

£uP] 



/ Tfl ' 

&A X ) = \/— TT ex P (*P r + «£i(p)*)v CT , (7) 

where £»(p) = y/mf + p 2 ; u„ and v a are orthonormal spin functions. 

The commutation rules for the operators aj jP and 6j iP are the same as for fermions. Henceforth, we neglect the spin 
interaction, and the spin parts in (7) are omitted. 

We define the vacuum state |0 > by: 

^ (i) (a;)|0>=0 (8) 

then 

< °Kp a tpl° >= !. 

< 0|6+ p ^ p |0 >= 1. 

The latter can be treated as a filled negative-energy sea. 

Substituting (5)-(7) into (4) and taking into account (8), we find: 

G°(n - r' l5 r 2 - v' 2 -E) = ^L 2 (p,q;£) exp(ip( ri - r[) + *q(r 2 - r 2 )), (9) 



where 



£i(p)£2 q 



£ - £i - £ 2 £ - £i + £ 2 £ + £i - £ 2 -E + £i + £ 2 

toito 2 AE(E 2 - £ 2 (p) - £ 2 (q)) 

ei(p)e 2 (q) S 4 - 2£ 2 (£ 2 + el) + (e 2 - e 2 ) 2 ' 



(10) 



Note that all the pole functions in G° (9) have equal positive-definite amplitudes. 

Thus, the equation (3) with (9) and (10) is the integral two-fermion wave equation. In the case when no constraint 
between ri and r 2 exists, this integral wave equation can be reduced to the following differential equation: 

£i(p)£ 2 (q) [E 4 - 2E 2 ( £l (pf + £ 2 (q) 2 ) + 
(£ 1 (p) 2 -£ 2 (q) 2 ) 2 ]e(r 1 ,r 2 ) = 

4mim 2 £ [E 2 - £l (p) 2 - £ 2 (q) 2 ] V{v x - r 2 )£(n,r 2 ), 

where p = — iV ri and q = — iV r2 . In the limit V — ■* 0, this differential equation has the four solutions E = 
±£i(p) ± £ 2 (q) for the plane wave functions (7). Here p and p are independent. These solutions correspond to the 
all poles of the propagator (9)-(10). So, one can see that the negative-energy part is not removed in the approach. 

To find solutions of the equation for the bound states, it is convenient to transform to coordinates in which the 
center-of-mass coordinate is separated out. 



III. THE EQUATION IN THE LABORATORY SYSTEM. 

The wave function of bound states can be presented in the form: 

0(n, r 2 ; E(g)) = £(r, g; E(g)) cxpfeR) (11) 

where r = ri — r 2 is the radius-vector between the particles; R is the radius-vector of the ccntcr-of-mass system and 
g is the wave vector of free motion of the bound system. 
The Fourier transform of the wave function is given by: 



V>(k, f;E)= J dridr 2 V(ri, r 2 ; E) exp(-ikn + tfr 2 ) (12) 



In (12) we can introduce new integration variables for fixed k and f: 



ri£i(k +r 2 e 2 f) , 1Q , 

r = ri-r 2 R = — — — — . (13) 

£l(k) +£ 2 (f) 



Substituting (13) into (12) and using (11), we obtain: 

^(k,f ;£(»))= fdrdB4(r, S ;E(g)) 



, .,. , ,-p . k£ 2 (f)+f £l (k) 

cxp(-,(k-f- g )R-,r ei(k)+g2(f) 



^3 JC , 1 , e „„/k£ 2 (f)+f£ 1 (k) 



(^(k-f-g)^ -^^^ ,,;^. (14) 

Now we calculate the (k, f) expansion term of Eq. (3). First, using the change of the variables (13) at the fixed p 
and q, we calculate the integral over r[ and r 2 in the right-hand side of Eq. (3). Then integrating with respect to p 
and q, we have: 

V(k, f ; E) = (27r) 3 (5(k - f - g)L 2 (k, f ; E) J dv' cxp {-{ S r')V{r')(;{r' , g), (15) 

where: 

From Eqs. (14) and (15) we obtain: 

£(s, g) = L 2 (f + g, f ; E{g)) J dv 1 cxp {-isv')V{r')Z(v\ g) (17) 

Given the solution of Eq. (16) in the form f (s, g), then from Eq. (17) we obtain the integral equation for the bound 

states: 



£(r, g) = J 7^3 J dr'L 2 (f(s, g) + g, f (s, S );E(g)) 



(27T) 

cxp(zs(r-r'))nr'K(r',g) (18) 

It can be easily shown that the integral equation (18) with account of (9) and (10) corresponds to the following 
differential equation for the eigenstates of the two-fermion system: 

£l (f (§, g) + g)£ 2 (f (§, g)) \E(gf - 2S(. 9 ) 2 (£ 1 (f (§, g) + g) 2 + £ 2 (f (s, g)) 2 )+ 



( £l (f(s,g)+g)^- £2 (f(s,g))^(r,g) = 

Am 2 E(g) [E(g) 2 - £l (f (§, g) + g) 2 - e 2 (f (S, g)) 2 )] U(r)£(r, g), (19) 

where f (s, g) is the solution of Eq. (16) and s = — iV r is the canonical momentum. 
One can check that substitution of the Fourier transform of the wave function: 



f ds' 

C(r,g) = / 7^T3-C(s', g) cxp (is'r) 



to Eq. (19) results in the integral equation (18). 

Note that Eq.(19) which can be used to study the two-fermion systems such as the electron-proton system and 
particle-antiparticle system, has a significantly non-Hamiltonian form. Unlike Hamiltonian equations, we cannot, 
in the general case, separate the contributions of the potential and kinetic energies to the binding energy of the 
eigenstates. The class of solutions for this equation can be distinguished from that for Hamiltonian equations. Below, 
we apply Eq. (19) to the electron-positron system. 

IV. THE EQUATION FOR THE ELECTRON-POSITRON BOUND EIGENSTATES IN THE 

CENTER-OF-MASS SYSTEM. 

For the electron-positron system mi = m 2 = m and q\ = — q 2 = e. We have g = in the center-of-mass system. 
Then from Eq. (16) we obtain that f = s, and Eq. (19) is reduced to the form: 

e(s) 2 E(0) \E 2 - 4e(s) 2 j £(r) = 4m 2 \e 2 - 2e(s) 2 j U(r)£(r), (20) 

where e(s) = \Jm? + s 2 and § = — zV r . 

Further we use the Coulomb interaction potential V(r) = — — (sec [1] and references therein). It is convenient 
to introduce dimcnsionless variables in Eq. (20): E — ► Em and r — ► A e x, where A e is the Compton length of the 
electron. Then Eq. (20) is reduced to the form: 

E{1 - A x )(£ 2 - 4 + 4A x )£(x) = -Aa(E 2 - 2 + 2A x )-£(x) (21) 

x 

First of all, we show that Eq. (20) contains the class of the positronium eigenstates of the Schrodinger equation 
with the Coulomb interaction potential, and with certain relativistic corrections. Using the replacement E = 2m + Eb, 
from (20) we obtain the effective Hamiltonian: 



s 2 

— + U(r) - 
m 


3_El 
4m 


E b s 2 t 2E b V(r) t s 4 
2m 2 m m 3 


s 2 U(r) 

2 

rrr 



H=- + V(v) -^L-^ + -^L^. + ^, (22) 

m Am lm z 

where we omitted five terms of the order of 0(a 6 ) and 0(a 8 ). 

One can see that Eq. (22) has the correct "unperturbed" Hamiltonian which includes the first and second terms 
in the right-hand side of (22). But the eigenvalues of (22) cannot agree to the order of ~ a 4 with the positronium 
spectrum [22] because we neglected the relatively small spin interactions. It is not essential for us here. Our aim is to 
show that Eq. (21) obtained with the use the free two-fermion propagator which satisfies the Lemann representation, 
has a new class of the electron-positron bound states. 

We find the s— symmetrical bound eigenstates in the form £(x) = exp (— j3x). Then from (21) we obtain: 

E{1 - [3 2 )[E 2 - 4(1 - P 2 )] exp (-fix) - 327r(a + f3E)5(x) + 

(2(3E[E 2 - 8(1 - f3 2 )] + Aa[E 2 - 2(1 - (3 2 )]) - exp (-/3x) = (23) 

For the positronium states, the delta-function term in Eq. (23) can be considered as a perturbation defined at x = 0. 
Omitting the perturbation (or considering Eq. (23) at x ^ 0), we find two algebraic equations in the unknowns E 
and /3: 



£(l-/? 2 )[£ 2 -4(l-/3 2 )]=0 (24) 

2(3E[E 2 - 8(1 - /3 2 )] + Aa[E 2 - 2(1 - f3 2 )} = 0. (25) 



Solution (24) yields: E\ = 2-^/1 — /3 2 , ft = 1 and £3 = 0. The latter, £3, does not correspond to a eigenstate 
because of divergence of L 2 (10) at E\ — e 2 in the center-of-mass system. 

Substituting E\ and ft into (25), we find the two eigenstates with the s— symmetry. This situation is only due to 
the non-Hamiltonian form of (20). The first eigenstate is the Is state of positronium: 



E 1 = 2* (l + y/1 - a 2 ) , ft = 2"*(l - y/l - a 2 ) . (26) 

In usual units the binding energy E b \ and exponential damping radius a,\ of the eigenstate are given by: 

£"61 = mc (2 - £1), 01 = a B — . 

Pi 

Here a^ is the Bohr radius. 

The second solution corresponds to the deep bound state: 

E 2 = -2a, & = 1. (27) 

In usual units the binding energy E b2 and exponential damping radius a 2 of this eigenstate are given by: 

E b2 = 2(1 + a)mc 2 , a 2 = A e . 

Now we take into account the delta-like term in Eq. (23) and solve this equation by the variational procedure. 
Multiplying (23) through by £(x) and then integrating with respect to x, we get the cubic equation: 



,3 , p2 4a/3 4(1 + 2/3 2 + 4/3 4 ) 8a/?(l + 3 ft 2 ) 

1 + ft l + f3 2 1 + ft 






solution of which gives us the three real roots Ei(/3), i = 1,2,3. The minima of the Ei(f3) functions were found by 
using Mathematica 2.2.1. 

The plot of Ei(/3) is shown on Fig. 1. One can see that the dependence for the first root has the minimum at 
Pi = 3.64769 x 10~ 3 . The corresponding value of the dimensionless total energy £i(/3i) = 1.99998669. Thus, the first 
solution is the Is state of positronium. In usual units we have the exponential damping radius of the wave function: 

a 
a\ = a B — = 2.00054 x a B 
Pi 

and the binding energy: 

E bl = mc 2 * (2 - £i(/3i)) = 6.80157eF. 

These values are in good agreement with the results (26). Thus, the delta-like term in Eq. (23) can really be considered 
as the small perturbation which results in the correction of the order of ~ a 4 to the positronium spectrum. 

The plot of •E l 2(p > ) is shown on Fig. 2. One can see that the dependence for the second root has the minimum at 
p2 = 0.725625. The corresponding value of the total energy ^(ft) = —7.94318 x 10~ 3 . This second solution is the 
deep bound state mentioned above. In usual units this eigenstate is determined by the damping radius: 

a 2 = ^ = 1.37812 x A e 

P2 

and the binding energy: 

E b2 = m*{2- E 2 ((3 2 )) = 2.00794 x mc 2 

Thus, the result obtained by the variational method, shows also the possibility of the unusual deep bound state for 
the electron-positron system. Comparing these values a 2 and E 2 with the result (27), we note that the delta- like term 
in Eq. (23) cannot be considered as the small perturbation for the deep bound state. Because the interaction between 
the particles at r — > is important, this state can, of course, be affected by the electron vacuum polarization. 
The third root E%{f}) of Eq. (28) has a wide maximum near (3 ~ and for this reason, it is not of interest. 



V. THE MASS OF THE BOUND SYSTEM FOR THE TWO TYPES OF THE EIGENSTATES. 

Mow we calculate the mass of the electron-positron system for both the positronium Is state and deep bound state 
£2- At low kinetic energies g <C m the solution of Eq. (16) is: 



Using (29), we have: 



f(s,g) = s-i g +^ + 0( 3 3 ). 



£ (f(s,g)+g) £ (f(s,g))= £ (s) 2 + - 9 2 ; 



£ (f (s, g) + g) 2 + £ (f (s, g)) 2 = 2 £ (s) 2 + \g 2 + MJ; 

2 e(s)^ 



(29) 



£ (f(s,g)+g) 2 - £ (f(s,g)) 2 = 2(sg). 
For E(g) we use the following expansion: 



(30) 



E(g) = E(0) 



2m s ' 



(31) 



where m s is the mass of the bound system at low kinetic energies; E(0) is the obtained total energy of the bound 
system at rest. 

Eqs. (29)-(31) are substituted in Eq. (19). Instead of f(r,g), we use the wave function £(r, 0) obtained for the 
system at rest. After simple manipulations and reduction of the equation by g 2 we get: 



<ei 



Ee(s) 2 (E 2 -2e(s) 2 ) + am 3 [3E 2 - 2e(s) 



Xe 



r 



\t>- 



2 < $\s 2 z (e 2 2 £ (§) 2 ) |C > -\E^ < Z\£ > +2E 2 < e| £ (s) 2 |C > + 



2am 3 E < f| — 1£ > +2am 3 E < £|-^-^|£ > . 

r ' £ (s) z r 



(32) 



Here s = — iV x and g = (0, 0, g) is proposed without loss of generality. 

The last term in the right-hand side of Eq. (32) can be transformed by using Eq. (20) which is written in the form: 



E < f I [E 2 - 4 £ (s) 2 ) |f >= -4am 3 < f | 



E 2 \ A e , 

e(s) z ) r 



(33) 



Using the above dimensionless values and taking into account Eq. (33), Eq. (32) is reduced to the form: 

2m 



<i\ 



\e 3 {\ - A x ) - E{\ - A x ) 2 + 2aE 2 - 
4 x 



\z>- 



4 7 . 1 4 



., 3 £ 2 - -E*) < ft > +-(1 - E 2 ) < e|A x |e > 



l<^^> + 1 -(l0aE-^)<^> 



(34) 



Finally, substituting f(x) = cxp (— fix) in Eq. (34), we get: 



_ 2m f£ 3 (l+/? 2 ) + 2£ 2 Q; 3-£(l + 2/? 2 + 5/3 4 ) 

m " W f (l0«£-8§)- ±£ 4 + ±(7 + 8/? 2 )-§(l + 2/? 2 )-f/3 4 

Substituting E\ and /3i which were obtained by the variational method, into Eq. (35), we arrive at: 

m s ,i = 1.99995121 x m (36) 

for the Is state of positronium. This value is in agreement with the positronium mass. 

The deep bound state is determined by E 2 — 7.94318 x 10~ 3 and j3 2 = 0.725625. Using these values, we get from 
Eq. (35): 

m s2 = -0.0097293 x2ra~ -9943.4eU. (37) 

Thus, the mass of the electron-positron system in the deep bound state is predicted to be very small in comparison 
with the positronium mass. 

VI. CONCLUSION 

The free single-fermion propagator used in QFT, is commonly adopted to be identical with the Green function for 
the Dirac differential equation. Using (1) and (2), this propagator is given by: 

G°(p;£0- 1P + m( l 1 



2e(p) \E-e(p) E + e(p) 

One can see the minus between the two pole functions. However, this propagator structure is common for Bose 
particles and does not satisfy the Lemann spectral representation for fermions [19,20]. 

The bare single-fermion propagator which satisfies the Lemann representation, can be written as [23] (the spin 
index is omitted): 

G°( P ; E) = -^ ( ^ + L^ (38) 

VP ' ; e(p) ^-e(p) E + e(p)J K ' 

with the plus between the two pole functions. 

This sign difference gives rise to the class of deep bound eigenstates (of course, in addition to the Bohr states) for 
the one-fermion problem in the Coulomb potential V(r) = — — , as shown in [23]. Using (38), we can show that the 
integral equation for the bound states: 

(f>(r;E) = f dr'G°(r-r';E)V(r')(f>(r';E) 

corresponds to the differential equation: 

e(p) \E 2 - e(p) 2 j 0(r; E) = 2mEV(r)<j)(r; E), 

which has the non-Hamiltonian form also. 

The use of Dirac's free propagator with the minus between the two pole functions causes the loss of the deep bound 
state solution for the one-fermion problem in the Coulomb potential [23] . This solution is lost because of the two pole 
functions are of the same order of magnitude for the deep bound state, unlike Bohr states. We have shown here that 
the same situation takes place for the electron-positron system as well. 

2 -. 

In the above consideration, the coupling constant is given by the fine structure constant a = f- = jkf. Thus the 
origin of these deep bound states distinguishes with the bound-state solutions of the two-fermion Dirac-like equations 
at the increasing of the coupling constant [1,7]. 

The wave equation (20) for the electron-positron system has non-Hamiltonian form. The two s— symmetrical 
solutions obtained exist only by this reason. But the eigenvalues for the positronium spectrum cannot agree to the 
order of ~ a 4 with well-known results [22] because we neglected the relatively small spin interactions. It was not 
essential for us here. We have shown that the use of the free two-fermion propagator which satisfies the Lemann 
representation, gives rise to the new class of the electron-positron bound states. The second solution obtained is the 
unusually deep bound state with the anomalous small mass of the bound system. Of course, this state can be affected 
by the electron vacuum polarization and spin interactions. 

Similar considerations can be carried out for the proton-antiproton system. We believe that this deep bound state 
with the anomalous small mass can also arise for the system. But one should take into account that these particles 
are not point ones, and vacuum polarization effect can be an essential factor. 
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FIG. 1. The /3 dependence of the first root Ex of Eq. (28). 
FIG. 2. The /3 dependence of the second root E 2 of Eq. (28). 



in 

O 



X 
i 



1.0 



0.5 



0.0 



-0.5 



-1.0 



-1.5 




0.0 0.002 0.004 0.006 0.008 



Fig. 1. 



A.I. Agafonov and E.A. Manykin. 

"Relativistic bound-state equations for two-fermion systems' 



-7.0 



-7.25 



en 
O 



x -7.5 

r N 



-7.75 



-8.0 




0.5 0.625 0.75 0.875 



1.0 



Fig. 2. 



A.I. Agafonov and E.A. Manykin. 

"Relativistic bound-state equations for two-fermion systems' 



